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ABSTRACT 



Synchronous binary asteroids may exist in a long-term stable equilibrium, where 
the opposing torques from mutual body tides and the binary YORP (BYORP) effect 
cancel. Interior of this equilibrium, mutual body tides are stronger than the BYORP 
effect and the mutual orbit semi-major axis expands to the equilibrium; outside of the 
equilibrium, the BYORP effect dominates the evolution and the system semi-major 
axis will contract to the equilibrium. If the observed population of small (0.1 - 10 km 
diameter) synchronous binaries are in static configurations that are no longer evolv- 
ing, then this would be confirmed by a null result in the observational tests for the 
BYORP effect. The confirmed existence of this equilibrium combined with a shape 
model of the secondary of the system enables the direct study of asteroid geophysics 
through the tidal theory. The observed synchronous asteroid population cannot exist 
in this equilibrium if described by the canonica l "monolithic" geophysic al model. The 
"rubble pile" geophysical model proposed by iGoldreich & Saril (|2009|) is sufficient, 
however it predicts a tidal Love number directly proportional to the radius of the aster- 
oid, while the best fit to the data predicts a tidal Lo ve number inverse l y prop ortional 



Goldreich & Sari (2009) models 



to the radius. This deviation from the canonical and 
motivates future study of asteroid geophysics. Ongoing BYORP detection campaigns 
will determine whether these systems are in an equilibrium, and future determination 
of secondary shapes will allow direct determination of asteroid geophysical parame- 
ters. 



Subject headings: Celestial mechanics — Minor planets, asteroids: general — Planets and 
satellites: dynamical evolution and stability — Planets and satellites: fundamental parameters 
— Planets and satellites: interiors 
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1. Introduction 



Synchronous binary asteroid systems may now provide a method of directly probing an 
asteroid's internal structure. They may be inhabiting a long-term stable equilibrium created by the 
opposing torques from mutual body tides and the binary YORP (BYORP) effect. From the tidal 
theory, this equilibrium would allow direct study of an ast eroid's geophysics for the first time. 



Goldreich & Saril (120091) . but not fully 



This potential equilibrium was suggested as a possibility by 
explored. 

For this long-term stable equilibrium to exist, the internal structure of binary asteroids 
must be different than th e usually assumed "monolithic" canonical internal structu re developed 



for astrophysical bodies (|Darwin 



1879 : 



Love 



1927b 



Goldreich 



1963 



Kaula 



19641) . Since the 



canonical theory was determined by studying much larger, hydrostatically balanced bodies 
under simplifying continuum mechanical approximations, it is not s urprising that this t heory 



Goldreich & Sari 



(2009) 



does not encompass small, very low gravity, "rubble pile" asteroids, 
derived an alternative theory for the internal structure by modifying the continuum mechanical 
approximations to include the effect of voids in the "rubble pile" structure. As they noted, this 
new internal structure does allow for this equilibrium to be stable in eccentricity. 

We study the internal structure of synchronous binary asteroid systems by making two 
observationally and theoretically motivated assumption s. First, these binaries formed via 



rotational fission and so share a common parent body ([ Margot et al 



Walsh et al. 


2008; 


Pravec et al. 


2010; 



Jacobson & Scheeres 



2011 



2002; 



Scheeres 



2007 : 



). Each component of the binary 



will share the same intensive properties such as density p, tidal dissipation number Q, rigidity fi, 
and yield strain ey. Extensive properties such as radius R, mass ratioj q, and tidal Love number k 



'Defined as the mass of the smaller (secondary) component divided by the mass of the larger 
(primary) component of the binary. 
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depend on the absolute sizes of the binary members. 



Second, these systems are observed proportional to the total time they spend in each state of 



their evolut ion. Binary asteroid 



popula tion (pvlargot et al. 



2002; 



systems comprise a sign ificant fracti on ( 



Pravec et al. 



2006). Both 



Cuk (2007) and 



5 ± 4%) of the NEA 



McMahon & Scheeres 



(|2010|) determined that if the BYORP effect dominates the evolution, then synchronous binaries 
can disrupt in much less than a million years, requiring frequent binary creation to maintain the 
observed population. If binaries are trapped in a long-term stable equilibrium, then binary creation 
could be infrequent. This conclusion is corroborated by evidence that the binary formation process 
from rotational fission is inefficient, requiring many rotational fission events per asteroid (taking 
many Y ORP timescales), in order to create a stable binary system that does not immediately 



disrupt (|Jacobson & Scheeres 



201 lb . The observed synchronous binary population may be 



residing in this equilibrium. 

After the formation of a stable binary asteroid system via rotational fission (or some other 
mechanism), tides immediately begin to dissipate energy from the system through the mechanical 
stressing of each body. The fastest process is the synchr onization of the secondary which occurs 



because of tides raised on the secondary by the primary (|Goldreich 



1963; 



Goldreich & Sari 



2009). 



This process is often referred to as tidal locking of the satellite. After the secondary has been 
synchronized, the system will evolve due to both tides and the BYORP effect, which requires a 
synchronous secondary. 



1.1. Tidal Evolution 

Relative motion between components in a binary system leads to tidal dissipation of energy 
and the transfer of angular momentum between spin and orbit states. Assuming spherical, 
homogenous bodies with identical compositions and a mutual orbit with low eccentricity, 



-5- 



an asteroid's first order geophysics can be characterized by two parameters: the tidal Love 
number and the tidal dissipation number. The tidal Love number k is the ratio of the additional 
gravitational potential produced by the redistribution of mass relative to the deforming potential, 
and theoretically is thought to depend on the size and internal properties of the body. The tidal 
dissipation number Q describes how effective the body is at tidally dissipating energy. It is a 
quality factor defined as: Q = 2nE / § ^ where E is the maximum energy stored in the 
tidal distortion during a cycle and § ^ dt is the energy dissipated over one cycle. 

As stated before, tides raised on the primary by the secondary synchronize the spin of the 



secondary to the mutua l orbit period; this is the fastest tidal evolutionary process (|Goldreich 



1963 



Goldreich & Sari 



2009J). Tidal evolution continues after synchronization of the secondary. As 
tides dissipate energy from the rotation states of binary members, the semi-major axis and the 
eccentricity of the mutual orbit change over time. Only after the rotation period of the primary is 
synchronous with the mutual orbital period as well and the system is thus doubly (completely) 
synchronous will binary systems not experience any tidal evolution. 

When the system is singly synchronous, tides raised on the primary by the secondary cause 
the semi-major axis to expand. From the first order theory, the semi-major axis time rate of change 
due to tides is: 

where a is the semi-major axis measured in primary radii R p , k p is the tidal Love number of 
the primary, co d = (At vGp/3) 1 ^ 2 is the surface d isruption spin limit for a sphere, and G is the 



gravitational constant (|Murray & Dermott 



1999). Due to the fast rotation of the primary compared 



to the mean motion, tidal evolution always expands the semi-major axis of the binary system. 

Tides raised on the primary by the secondary cause the eccentricity to grow, while tides 
raised on the secondary by the primary cause eccentricity to damp. From the first order theory, the 
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tidal eccentricity time rate of change due to tides is: 

. _ 57k p q 1/3 - 84k s / bj d e \ 
6T ~ 8Q WW 



2/3, 



Q 



(2) 



where e is the eccentricity and k s is the tidal Love number of the secondary ([Murray & Dermott 
19991) . Tidal evolution can cause excitation or damping of the eccentricity depending on the 
system mass ratio and the tidal Love numbers of each component. 

There are two developed theories for the internal structure of asteroids: "monolith" and 
"rubble pile." Evidence that asteroids have a "rubble pile" internal structure rather than a 
"monolithic" interior includes: the Hayabusa mission to Itokawa showing no obvious impact 
craters and the appearance o f a structure made entirely from shattered frag ments of different siz e 



scales (|Fujiwara et al 



2006|) . numerical modeling of collisions of asteroids (|Mic 



lei et al 



mass and volume measurements from the NEAR Shoemaker 



1997|) and radar observations of 1999 KW 4 (|Ostro et al 



lyby of Mathilde (lYeomans et al 



2001), 



2006) showing mean densities that are 



lower than their constitutive elements, the rotational speed limit period (~ 2.2 hours) observed 
amongst asteroids with diameters larger than ~ 2 00 m, which corresp onds to the critical disruption 



spin rate of a self-gravitating, strengthless body (IPravec et al. 



20071), and that rotationa 



of strengthless bodies is responsible for the asteroid pair population (jPravec et al 



fission 



20101). All of 



this evidence suggests that a "rubble pile" internal structure is a more realistic assessment than a 
"monolith" internal structure, however for completeness we will consider both theories. The two 
theories become distinct when assessing the functional form of the tidal Love number. 

The dimensionless rigidity of the body fl, which can be thought of as the ratio of the fluid 



strain to the elastic strain, can be used to defined the tidal Love number: k 



this a pproximation is true when the fluid strain dominates the elastic strain (jGoldreich & Sari 



.5/(1 + u) w 1.5/ /2, 



2009b . According to the canonical "monolith" theory the dimensionless rigidity has the form 
fl = 19 fi/(2ui^pR 2 ), where fx is the rigidity of the body. To first order, the tidal Love number for a 
"monolith" depends on the size of the body as k M oc R 2 . 
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Goldreich & Sari 



(|2009l) developed an alternative "rubble pile" tidal Love number theory 
by studying how introducing voids increases the stress across contact areas. The dimensionless 
rigidity of a "rubble pile" is smaller t han that of a "monolith" of the same size. Using 



conceptual and dimensional arguments, 



Goldreich & Sari 



(|2009r) determine that the "rubble pile" 



dimensionless rigidity should scale as Jir > (/i/ey) 1 / 2 where ey is the yield strain. To first order, 
the tidal Love number for a "rubble pile" depends on the size of the body as k R oc R. Both of 
these theories will be explored in later sections. 



1.2. BYORP Evolution 



Th e BYORP effect is t he summation of radia t ive ef fects on a synchronous sec 



ondary (|Cuk & Burns 



2005). 



McMahon & Scheeres 



3) 



(120101) using averaging theory showed that 



this effect secularly evolves both the semi-major axis and the eccentricity, because radiative forces 
acting on asymmetries in the shape of the secondary create torques on the mutual orbit that persist 
after averaging over the mutual orbit, the heliocentric orbit, and the precession of the node. To 
first order, the evolution of the semi-major axis and the eccentricity only depends upon a single 
constant term B t hat represents the averaged a cceleration in the direction parallel to the motion 



of the secondary (|McMahon & Scheeres 



1201(1) . The BYORP coefficient B depends only on the 
shape of the secondary, and can be thought of as a ratio relating the asymmetric area, on which 
the average force of the BYORP effect acts, to the total area. A symmetric body-a sphere or an 
ellipsoid-has a value: B = 0. The BYORP coefficient has a maximum magnitude: B = 2 in either 
direction, but commonly has a small value: B ~ 10~ 3 . The BYORP effect can either expand or 
shrink the semi-major axis, with the sign of the eccentricity evolution always opposite that of the 
semi-major axis evolution. From the first order theory, the evolution of the semi-major axis and 
eccentricity is: 



a B 



2vr 



\u d pR$J 



Q 



1/3 
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=F 



3H Q B ( a^e \ JTT~g 



^7T 



\UdpRj) 



1/3 



(4) 



where H & = F /(a|^/l — e|), F Q is the s olar radiation constant, and and e are the 



heliocentric semi-major axis and eccentricity (jMcMahon & Scheeres 



2010). 



2. Joint Evolution 

Synchronous binary asteroids will evolve under the influence of both tides and BYORP. 
There are two scenarios: joint expansive evolution and joint opposing evolution depending on 
the direction of the BYORP torque, which nominally has an equal chance of pointing in either 
direction. Both tides and the BYORP effect can change the energy of the system over time, and 
the BYORP effect can also changes the system angular momentum. It is convenient to define a 
parameter A, which compares the strengths of these two effects on the semi-major axis of the 
system, and E which does the same for the eccentricity: 

A = JM = 2?r ^p fc p^g 4/3 (5) 

| Ob | BH Q Qa 7 

E= \p\ji»-nK\ A 

\ e B\ 2 

where K = k s / (kpq 1 ' 3 ) is an important tidal Love number relation. For the "monolith" model 
k oc R 2 , so K = g 1 / 3 , and for the "rubble pile" model k oc R, so K = 1. If A < 1 then the 
BYORP effect dominates semi-major axis evolution, A > 1 then tides dominate the semi-major 
axis evolution, and if A = 1 then the two are balanced and the system semi-major axis will not 
evolve via either effect. Which effect dominates eccentricity is shown similarly by E. 



2.1. Joint Expansive Evolution 



During joint expansive evolution, both effects are growing the semi-major axis. Tides are 
removing energy by spinning down the primary and transferring angular momentum from the 
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primary spin to the orbit, while the BYORP effect is directly adding both energy and angular 
momentum to the orbit. The eccentricity of the mutual orbit can be excited or damped depending 
on the system mass ratio and the tidal Love numbers of the components. The condition for 
stability or damping in eccentricity is: 

19 - 28K A 

A < 1 (7) 

For "monoliths," K = q 1 ^ 3 and so low mass ratio systems can have growing eccentricity, if 
q < (19/28) 3 and A > 2/(19 - 28g 1/3 ). For "rubble pile" systems, K = 1 and thus the condition 
is always satisfied, and so the eccentricity of binary asteroids with "rubble pile" internal structures 
will always damp. 

Since the semi-major axis growth is always positive, systems that damp eccentricity will 



gr ow towards the Hi l 



of 



Cuk & Nesvornv 



sphere and eventually disrupt, disregarding the possible evolutionary paths 



(|2010l) . Systems that grow in both semi-major axis and eccentricity will do 
so until the evolutionary equations in the first order theory no longer apply and more complicated 
dynamics become dominant. 



2.2. Joint Opposing Evolution 

During joint opposing evolution, tides are acting to grow the semi-major axis, but the BYORP 
effect is acting to shrink it. Thus the system evolves towards an equilibrium point, where these 
two effects balance. An equilibrium exists for the evolution of the semi-major axis at 

1/7 

(8) 



' 27rk p u 2 d pR 2 p q^ 



BH Q Q 

This is a stable equilibrium, and regardless of the initial semi-major axis, the system will evolve 
in semi-major axis to this equilibrium point. Tides are still transferring angular momentum to 
the orbit and removing energy from the system, but now the BYORP effect is removing angular 



momentum and energy from the system. Interior to the equilibrium semi-major axis, tides 
dominate and the orbit grows due to the increase of angular momentum to the orbit. Exterior to 
the equilibrium, the BYORP effect controls the evolution and the orbit shrinks due to the decrease 
of orbital angular momentum. At the equilibrium the amount of angular momentum removed 
from the orbit equals the amount tidally transferred into the orbit from the primary, de-spinning 
it. The torque on the primary is equal in strength to the BYORP torque T B and so the rate of the 
de-spinning of the primary is: 

Y B A5H Q B ( a p \l + q 



The eccentricity of the mutual orbit can be excited or damped depending on the system mass 
ratio and the tidal Love numbers of the components. The condition for stable evolution without 
growth in eccentricity is: 

±^*4<-l (10) 

For "monolith" asteroids, K = q 1 ^ 3 and so it is possible for any mass ratio system to be unstable in 
eccentricity if A < 2/(28g 1//3 — 19) and in fact, any system with a mass ratio q < (19/28) 3 « 0.31 
will be unstable regardless of A. At the equilibrium, A = 1 and so binaries with "monolith" 
internal structures will grow in eccentricity if the system mass ratio q < 27/64 « 0.42. This 
would exclude the observed synchronous binaries from existing in this equilibrium condition, 
since they all have mass ratios q < 0.2 as shown in Table [Q Therefore these systems might enter 
equilibrium but could not remain there for a long period of time since their eccentricities would 
continue to grow. 

For binary asteroids with "rubble pile" internal structures, K = 1 and the stability condition is 
always satisfied, and thus the eccentricity will always be damped. Thus the observed synchronous 
binary population can exist in the equilibrium without growth in eccentricity. 
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3. Implications for the Synchronous Binary Asteroid Population: 



If the observed synchronous population is assumed to be in this joint opposing evolutionary 
equilibrium state, then A = 1 and Eqn. [5] is solved for the three unknown quantities (B, Q, 
and k p ): 

BO 2ixuj 2 ,pRla^ 

(ID 



k r , 



Hc,a 7 



Table Q] lists and Fi g. [TJ plots BQ/ k p for each of the known s ynchronous binary systems using 



observational data (PrayecetaL 



2006; 



Pravec & Harris 



2007|) . The tidal dissipation number Q is 



an intensive property that we expect to be similar for a 



been estimated to be Q = 10 2 (|Goldreich & Soter 



1 of these bodies, and for small bodies has 



19661) . The scatter and size dependence in Fig.Q] 



should be from B and k p , respectively. The BYORP coefficient B does not depend on size; B 
does depend on the shape of the secondary. Asteroid sha pes can vary greatly introducing scatter 



in the BYORP coefficient. 



McMahon & Scheeresl (120 10h estimate B = 10 3 from the shape 



model of the secondary of 1999 KW 4 , which does not appear symmetric. The BYORP coefficient 
may vary over a few orders of magnitude, especially towards smaller values corresponding to 
sec ondaries that are more sy mmetric. The tidal Love number may have a dependance on size, 



and 



Goldreich & Sari 



(|2009l) predict that the tidal Love number k p = 10 5 R P for a "rubble pile" 



internal structure. 

In the top plot of Fig. [H the solid line plots a simple theoretical model of BQ/k p using the 
estim ates of each value f rom ab ove, so that BQ/k p = lO 4 /?" 1 . Fitting the proportionality constant 



of the 



Goldreich & Sari 



(|2009|) tidal Love number model does not significantly change the results. 
Dashed lines indicate scatter from the BYORP coefficient (one order of magnitude larger and two 
orders of magnitude smaller). This model works acceptably well for systems with primary radii 
R p = 2 km but predicts values too large for systems with much smaller primaries. Deviations 
away from the model appear to be correlated with primary size, and so this may indicate that the 
tidal Love number has the wrong primary radius dependance. 
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0.2 0.5 1.0 2.0 5.0 



Primary Radius [km] 



10 5 




0.2 0.5 1.0 2.0 5.0 



Primary Radius [km] 



Fig. 1. — BQ/k p were calculated directly from observed quantities according to equation ITTI for 
each known synchronous binary, and plotted as a function of primary radius R p along with 1-sigma 
uncertainties (for clarity, the same data are shown in both plots). This d ata is also listed in Table 



Goldreich & Sari 



(2009) 



1. The solid line in the top plot shows the tidal Love number model from 
for asteroids with "rubble pile" internal structures: k p = 10~ 5 R P . The bottom plot is a fitted model 
to the data: k p = 2.5 x 10" 5 /?" 1 . For both models, the tidal dissipation number Q = 10 2 and the 
BYORP coefficient B = 10~ 3 is assumed. The dashed lines indicate the range of predicted scatter 
in the model due to the BYORP coefficient (possibly 10 times stronger or 100 times weaker). 
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If a power law is fit to the logarithmic data, then the best fit is BQ/k p = 4 x 10 3 i? p (shown 
as the solid line in the bottom plot of Fig.Q]). Using the same models for B and Q as above, then 
the tidal Love number dependence is k p = 2.5 x 10 -5 /?" 1 . Systems with this tidal Love number 
dependence are stable in eccentricity at the long-term equilibrium in semi-major axis. Deviations 
from this model do not appear to have a dependance on primary radius, and furthermore the 
scatter follows our expectations for scatter due to the BYORP coefficient, namely that the largest 
deviations are about two orders below, while most systems fall within an order below and above. 



4. Conclusions 

Small, synchronous binary asteroids may be residing in a long-term stable equilibrium. 
To inhabit the equilibria, asteroids cannot have a "monolith" internal structure but must have 
"rubble pile" interior. The tidal Love number is expected to vary with size, however the data 



Goldreich & Sari 



(2009) 



suggest k oc RT 1 , which is different than the kn oc R predicted by the 
"rubble pile" theory. This may highlight the difficulty of using a modified continuum theory to 
model "rubble piles." This suggests that the tidal Love number theory is incomplete, and future 
work should examine this closely. 

Current ca mpaigns to detect the BYORP effect may not achieve th e results predicted by 



previous studies (|Pravec & Scheirich 



2010; 



McMahon & Scheeres 



2010J). Tides may have stopped 



the synchronous population from evolving long ago. Existence of a long-lived equilibrium 
prevents rapid disruption due to the BYORP effect, and so binary formation at a furious rate is no 
longer required to maintain a steady-state synchronous binary population. It may take many years 
of observation in order to conclude that systems are in this equilibrium and are truly no longer 
evolving. 



In the future, if B is determined through secondary shape modeling for each system, the 
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geophysical parameters Q/k could be assessed direcdy. 
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Table 1. Properties and calculated BQ/k p of synchronous binary asteroid systems. 
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Note. — Data with 1-sigma uncertainties taken from the April 8, 2010 binary asteroid param- 



eter release from |http : / /www .asu.cas. cz/ -asteroid/binastdata .htm as compiled 



by methods and assumptions described in 



Pravecetal 



( 2006 ) and 



Pravec & Harris 



(2007). 



- 17- 

REFERENCES 

Cuk, M., & Burns, J. A. 2005, Icarus, 176, 418 

Cuk, M. 2007, ApJ, 659, L57 

Cuk, M., & Nesvorny, D. 2010, Icarus, 207, 732 

Darwin, G.H. 1879, Philos. Trans. R. Soc. London 170, 135 

Fujiwara, A., et al. 2006, Science, 312, 1330 

Goldreich, P. 1963, MNRAS, 126, 257 

Goldreich, P., & Soter, S. 1966, Icarus, 5, 375 

Goldreich, P., & Sari, R. 2009, ApJ, 691, 54 

Jacobson, S. A., & Scheeres, D. J. 2011, Icarus, DOI: 10. 1016/j.icarus.201 1.04.009 

Kaula, W. M. 1964, Reviews of Geophysics and Space Physics, 2, 661 

Love, A.E.H. 1927, A Treatise on the Mathematical Theory of Elasticity. Dover, New York 

Margot, J. L., Nolan, M. C., Benner, L. A. M., Ostro, S. J., Jurgens, R. R, Giorgini, J. D., Slade, 
M. A., & Campbell, D. B. 2002, Science, 296, 1445 

McMahon, J., & Scheeres, D. 2010, Celestial Mechanics and Dynamical Astronomy, 106, 261 

McMahon, J., & Scheeres, D. 2010, Icarus, 209, 494 

McMahon, J., & Scheeres, D. 2010, Bulletin of the American Astronomical Society, 42, 1093 
Michel, P., Benz, W., Tanga, P., & Richardson, D. C. 2001, Science, 294, 1696 
Murray, C. D., & Dermott, S. F. 1999, Cambridge University Press 



- 18- 

Ostro, S. J., et al. 2006, Science, 314, 1276 
Pravec, P., et al. 2006, Icarus, 181, 63 

Pravec, P., Harris, A. W., & Warner, B. D. 2007, IAU Symposium, 236, 167 
Pravec, P., & Harris, A. W. 2007, Icarus, 190, 
Pravec, P., et al. 2010, Nature, 466, 1085 

Pravec, P., & Scheirich, P. 2010, Bulletin of the American Astronomical Society, 42, 1055 
Scheeres, D. J. 2007, Icarus, 189, 370 

Walsh, K. J., Richardson, D. C, & Michel, P. 2008, Nature, 454, 188 
Yeomans, D. K., et al. 1997, Science, 278, 2106 



This manuscript was prepared with the AAS IATgX macros v5.2. 



